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ABSTRACT 

The Stokes parameters of the pulsed synchrotron radiation produced in the 
striped pulsar wind model are computed and compared with optical observations 
of the Crab pulsar. We assume the main contribution to the wind emissivity 
comes from a thin transition layer where the dominant toroidal magnetic field 
reverses its polarity. The radial component of the field is neglected, but a small 
meridional component is added. The resulting radiation is linearly polarized 
(Stokes V = 0). In the off-pulse region, the electric vector lies in the direction 
of the projection on the sky of the rotation axis of the pulsar. This property is 
unique to the wind model and in good agreement with the data. Other prop- 
erties such as a reduced degree of polarization and a characteristic sweep of the 
polarization angle within the pulses are also reproduced. These properties are 
qualitatively unaffected by variations of the wind Lorentz factor, the electron 
injection power law index and the inclination of the line of sight. 

Subject headings: MHD — Plasmas — Polarization — Pulsars: general — Radi- 
ation mechanisms: non-thermal 



1. Introduction 

The high-energy, pulsed emission from rotating magnetized neutron stars is usually 
explained in the framework of either the polar cap (Sturrock 1971; Ruderman & Sutherland 
1975) or the outer gap models (Cheng et al. 1986; Hirotani & Shibata 1999), for a review 
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see (Harding 2001). Although the existence of such gaps is plausible (Petri et al. 2002), these 
models still suffer from the lack of a self-consistent solution for the pulsar magnetosphere 
and are based on the assumption that the magnetic field structure is that of a (Newtonian) 
rotating dipole. Nevertheless, recent observations of the polarization of the optical pulses 
from the Crab (Kellner 2002; Kanbach et al. 2003) motivated detailed comparative studies 
of the emission from these models (Dyks et al. 2004) as well as a new variant of the outer gap 
introduced by Dyks & Rudak (2003) and called the "two-pole caustic" model. This model 
has a gap extending from the light-cylinder to the polar caps. It reproduces the twin-peak 
pulse intensity profiles by associating each peak with a different magnetic pole. In all of 
these models, the radiation is produced within the light cylinder. However, the pulse profile 
is determined by the assumed geometry of the magnetic field and the location of the gaps. 
Despite earlier claims (Romani & Yadigaroglu 1995), it now appears that neither the polar 
cap nor the outer gap model are able to fit the optical polarization properties of the Crab 
pulsar. The two-pole caustic model, on the other hand, provides a qualitatively reasonable 
fit (Dyks et al. 2004). 

An alternative site for the production of pulsed radiation has been investigated recently 
(Kirk et al. 2002), based on the idea of a striped pulsar wind, originally introduced by 
Coroniti (1990) and Michel (1994) and elaborated by Lyubarsky & Kirk (2001) and Kirk 
& Skjaeraasen (2003). Compared to the inner, outer and slot gaps, much less progress has 
been made in understanding the particle acceleration problem in this model. As proposed 
by Lyubarskii (1996), it is assumed that magnetic energy is released by reconnection in 
the thin regions where the toroidal field reverses its polarity, producing a non thermal elec- 
tron/positron population. But, unlike in Lyubarsky's model, emission from the striped wind 
originates outside the light cylinder and relativistic beaming effects are responsible for the 
phase coherence of the synchrotron radiation. A strength of this model is that the geometry 
of the magnetic field, which is the key property determining the polarization properties of 
the emission, is relatively well-known. Because relativistic, magnetically dominated winds 
seem to collimate only very weakly (Bogovalov 2001; Vlahakis 2004), the dominant field 
component well outside the light-cylinder is toroidal, irrespective of the geometry of the 
magnetic field near the stellar surface. 

In this Letter, we use an explicit asymptotic solution for the large-scale field structure 
related to the oblique split monopole and valid for the case of an ultrarelativistic plasma 
(Bogovalov 1999). This is combined with a crude model for the emissivity of the striped 
wind and of the magnetic field within the dissipating stripes themselves. We calculate the 
polarization properties of the high-energy pulsed emission by generalising a method applied 
to gamma-ray burst models by Lyutikov et al. (2003), and discuss and compare our results 
with optical observations of the Crab pulsar. 
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2. Calculation of the Stokes Parameters 

Our magnetic field model is based on the asymptotic solution of Bogovalov (1999), valid 
for r ^> t"l, and modified to take account of a finite width of the current sheet and with a small 
additional meridional component in the sheet, which serves to prevent the magnetic field from 
becoming identically zero. In order to simplify the Lorentz transformations, the small radial 
component of the field is neglected. In spherical polar coordinates (r, 9, ip) centered on the 
star and with axis along the rotation axis, the radial field is small (B r ~ B L r1/r 2 ) and the 
other components are: 

{B e , B^} 
r} e (A e ,r,9,ip,t) 

Here, B^ is a fiducial magnetic field strength, v is the (radial) speed of the wind, f2* = c/Vl 
is the angular velocity of the pulsar, with the radius of the light cylinder, a is the angle 
between the magnetic and rotation axes, bi t2 are parameters controlling the magnitude of the 
meridional field in the two current sheets present in one wavelength, and A e >ip are parameters 
quantifying the sheet thickness. The functional form of B^ is motivated by exact equilibria 
of the planar relativistic current sheet (see Kirk & Skjaeraasen 2003). However, in these 
equilibria the Bg component, which has an important influence on the polarization sweeps, 
is arbitrary. The Bg we adopt corresponds to a small circularly polarized component of the 
pulsar wind wave, such as is expected if the sheets are formed by the migration of particles 
within the wave, as described qualitatively by Michel (1971). 

For the particle distribution, we adopt an isotropic electron/positron distribution given 
by N(E,p,r,t) = K(r,t) E~ p where K(r,t) is related to the number density of emitting 
particles. The radial motion of the wind imposes an overall 1/r 2 dependence on this quantity, 
which is further modulated because the energization occurs primarily in the current sheet. 
The precise value in each sheet is chosen to fit the observed intensity of each sub-pulse. In 
addition, a small dc component is added, giving the off-pulse intensity. For the emissivity, we 
use the standard expressions for incoherent synchrotron radiation of ultrarelativistic particles 
in the Airy function approximation (Ginzburg & Syrovatskii 1969; Melrose 1971). Following 
Kirk et al. (2002), we assume the emission commences when the wind crosses the surface 
r = r > r L . 

The calculation of the Stokes parameters as measured in the observer frame involves 
simply integrating the emissivity over the wind. However, this requires special care, because 



B h — {b h2 Vei^e, r, 0, <p, t), rj v (A v , r, 0, <p, t)} 



tanh 



A v ^cos 9 cos a + sin 9 sin a cos j</? — f2* (t j j j 



1 d Vip (A e ,r,9, V ,t) 

Ag dif 



-4- 



the Lorentz boost from the rest frame of the emitting plasma involves not only beaming and 
Doppler shift, but also a change in the polarization angle due to the effects of aberration. 
Lyutikov et al. (2003) performed this calculation in the context of a gamma-ray burst model, 
assuming a relativistic shell of emitting plasma containing only a toroidal magnetic field. 
However, the relatively simple form of the emissivity function in that case means that two of 
the four Stokes parameters integrate to zero: U — V — 0, corresponding to linear polarization 
with constant position angle. Here we extend this method by adding a B e component and 
allowing for a more general, space and time dependent emissivity. 

After straightforward but lengthy manipulations involving Lorentz transformations, we 
find the Stokes parameters as measured by an observer at time t Q b s are given by the following 
integrals: 

P+7/3 



w r+oo pit pZtt p+ l 

Qlo ( (^obs) / / / s (r,9,<f,t Tet ) < cos (2%) > r 2 sm 9 dr d9 dip (2) 

U u J Jr ° h Jo { sin(2x) J 

where the retarded time is given by t iet = t Q b s + n ■ f/c and n is a unit vector along the 
line of sight from the pulsar to the observer. In this approximation the circular polarization 
vanishes: V = 0. The function s is defined by: 



P+3 



b r ——rA~ 



s (r,9,<P,t) = KK(f,t)u-— V— ^-y/l -(Vn-b) 2 j (3) 

where uo is the (angular) frequency of the emitted radiation, and k is a constant factor that 
depends only on the nature of the radiating particles (charge q and mass m) and the power 
law index p of their distribution: 

VSl f Sp + 7 \ { 3p-l \ \q\ 3 /3|g|\^ 
K ~ 2^4 rE A^~J rEu ^^2-J I^Tc {ri^J (4) 

— * 

with r Eu the Euler gamma function and V the Doppler boosting factor V = 1/T (1 — (3 ■ ft). 
The direction of the local magnetic field in the observer's frame is given by the unit vector 

— * 

b and the simplifying assumption has been made that this field has no component in the 
direction of the plasma velocity: b ■ f3 — 0. (In this case the magnetic field transformation 
from the rest frame B' to the observer frame B is just B' = B /Y and, thus, its direction 
remains unchanged.) The angle x measures the inclination of the local electric field with 
respect to the projection of the pulsar's rotation axis on the plane of the sky as seen in 
the observer's frame. The degree of linear polarization is defined by n = \JQ 2 + U 2 /I. 
The corresponding polarization angle, defined as the position angle between the electric field 
vector at the observer and the projection of the pulsar's rotation axis on the plane of the 
sky is x — 1/2 arctan(£7/Q). 
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3. Results 

Using a model similar to that described above, Kirk et al. (2002) computed the total 
radiation intensity and compared the spacing of the resulting twin-pulse profile with ob- 
servations of the Crab pulsar. They found an obliquity a = 60° and an inclination of the 
rotation axis to the line of sight £ = arccos(n • = 60°. In the following, we adopt 

these parameters, and set the radius at which emission switches on to be r = 30?"l. 

In the model of Kirk et al. (2002), the current sheet was assumed to be thin, which 
results in sharp profiles with very similar shapes for each of the subpulses. The upper left 
panel of Fig. 1 shows the intensity (Stokes parameter J) computed using our smoothed 
profile with A e — 1, A v = 5, b\ = 0.1 and b 2 = 0.08 for each subpulse. The electron 
density is K = [{(r L B L )/(r B)Y P+l ^ 2 + e - l]/[r 2 (1-0.6 77^] where the parameter e = 0.05 
sets the minimum electron density between the current sheets (in normalized units). The 
denominator (1 — 0.6779) introduces an asymmetry in the relative pulse peak intensity. The 
variation of the magnetic field and the particle density along the line of sight, are shown in 
the bottom panels of Fig. 1. 

The upper panels of Fig. 1 show the results of our computations on the left and the 
corresponding observed quantities (Kellner 2002; Kanbach et al. 2003) on the right. It should 
be noted that these data are preliminary. In particular, the measured degree of polarization 
may be subject to revision (Kanbach priv. comm.). Comparison with the upper right-hand 
panel shows that the model reproduces the observed pulse profile quite accurately. However, 
the idealised transverse geometry of the magnetic field in the (presumably turbulent) sheet 
leads to a rapid variation in phase of the term in parentheses on the right-hand side of 
Eq. (3), giving rise to a small notch-like feature visible in the peak of the sub-pulse, that 
proves difficult to eliminate. In this example, we adopted an electron power law index 
p = 2, as suggested by the relatively flat spectrum displayed by the pulsed emission between 
optical and gamma-ray frequencies (Shearer & Golden 2001; Kanbach 1998; Kuiper et al. 
2001). Results are shown for two values of the Lorentz factor of the wind: T = 20 (solid line) 
and r = 50 (dotted line). For convenience, the maximum intensity is normalized to unity. 
The timescale is expressed in terms of pulse phase, corresponding to the initial time t = 
and 1 to a full revolution of the neutron star and thus one period t = 27r/fi*. 

The degree of polarization is shown in the two middle panels of Fig. 1. According to our 
computations (left-hand panel) this displays a steady rise in the initial off-pulse phase, that 
steepens rapidly as the first pulse arrives. During the pulse phase itself, the polarization 
shrinks down to about 10%. Theoretically the maximum possible degree of polarization 
is closely related to the index p of the particle spectrum. In the most favorable case of a 
uniform magnetic field, it is given by II max = (p + l)/(p + 7/3). 
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However, in the curved magnetic field lines of the wind, contributions of electrons from 
different regions have different polarization angles. Consequently, they depolarize the overall 
result when superposed. We therefore expect a degree of polarization that is at most Il max . 
For the example shown in figure 1, II max (p = 2) = 69.2%, well above the computed value, 
which peaks at 52%. 

The lower panels of Fig. 1 show the polarization angle, measured against celestial North 
and increasing from North to East. Our model predicts this angle relative to the projection 
of the rotation axis of the neutron star on the sky, which we take to lie at a position angle 
of 124°, following the analysis of Ng & Romani (2004). In the off-pulse stage, the electric 
vector of our model predictions lies almost exactly in this direction, since it is fixed by the 
orientation of the dominant toroidal component B v of the magnetic field. Such a relation 
between the off-pulse angle of polarization and the toroidal magnetic field at large distance 
(close to the light cylinder) was already suggested by Smith et al. (1988). In the rising phase 
of the first pulse, B v decreases, whereas Bq increases, causing the polarization angle to rotate 
from its off pulse value by about 50°, for the chosen parameters. However, for a weaker Bq 
contribution, as in the second pulse, the swing decreases. This effect can also be caused by 
a relatively large beaming angle, (i.e., low Lorentz factor wind). The basic reason is that 
contributions from particles well away from the sheet center are then mixed into the pulse, 
partially canceling the contribution of the particles in the center of the current sheet, which 
favor x — 124° ± 90°, and enforcing x — 124°. On the other hand, a very high value of the 
Lorentz factor or large values of &i ;2 reduce the off-center contribution, leading, ultimately, 
to the maximum possible 90° sweep between off-pulse (-B^-dominated) and center-pulse (B e - 
dominated) polarizations, followed by another 90° sweep in the same sense when returning 
to the off-pulse. Thus, in general, in the middle of each pulse, the polarization angle is 
either nearly parallel to the projection of the rotation axis, or nearly perpendicular to it, 
depending on the strength of Bq and on T. This interpretation is confirmed by computations 
with T = 50 that show a larger sweep, as shown in Fig. 1. 

The optical polarization measurements suggest that in the centre of the pulses the 
position angle is close to 124°. In the declining phase of both pulses, the angle reaches a 
maximum before returning to the off-pulse orientation. Note that in both cases the swing 
starts in the same direction, (counterclockwise in figure 1). This is determined by the 
rotational behavior of the Bg component, implying that this changes sign between adjacent 
sheets, as in Eq. (1). The observed off-pulse position angle is closely aligned with the 
projection of the rotation axis of the pulsar, in accordance with the model predictions. 

In addition to models aimed at providing a framework for the interpretation of the 
emission of the Crab pulsar, we have performed several calculations with different Lorentz 
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factors T, injection spectrum of relativistic electrons p and inclinations of the line of sight £. 
The general characteristics of the results are: For low Lorentz factors, independent of p and £, 
the relativistic beaming becomes weaker and the pulsed emission is less pronounced, because 
the observer receives radiation from almost the entire wind. For instance, taking T = 2 and 
p = 2 or 3, the average degree of polarization does not exceed 20 % and the swing in the 
polarization angle is less than 30°. For high Lorentz factors T > 50, the strong beaming 
effect means that the observer sees only a small conical fraction of the wind. The width of 
the pulses is then closely related to the thickness of the transition layer. The degree of linear 
polarization flattens in the off-pulse emission while it shows a sharp increase followed by a 
steep decrease during the pulses. Due to the very strong beaming effect, only a tiny part of 
the wind directed along the line of sight will radiate towards the observer. In the off-pulse 
phase, the polarization angle is then dictated solely by the component "attached" to 
the line of sight, and the degree of linear polarization remains almost constant in time. For 
very high Lorentz factors, the behavior of polarization angle and degree remain similar to 
those of Fig. 1, with perfect alignment between polarization direction (electric vector) and 
the projection of the pulsar's rotation axis on the plane of the sky in the off-pulse phase 
and two consecutive polarization angle sweeps of 90° in the same sense during the off-pulse 
to center-pulse and center-pulse to off-pulse transitions. This mirrors the fact that emission 
comes only from a narrow cone about the line of sight of half opening angle 9 ps 1/T. 

For given values of T and £, the particle spectral index p affects only the average degree 
of polarization degree but not the light curve nor the polarization angle. For example, taking 
T = 10 and £ = 60°, a spectral index of p = 2 leads to an average polarization of LT = 19.2% 
whereas for p = 3 it leads to IT = 30.8%. 

4. Conclusions 

In the striped wind model, the high energy (infra-red to gamma-ray) emission of pulsars 
arises from outside the light cylinder, in accordance with the early suggestions of Pacini & 
Rees (1970) and Shklovsky (1970). It provides an alternative to the more intensively studied 
polar cap, outer gap and two-pole caustic models. All models contain essentially arbitrary 
assumptions concerning the configuration of the emission region and the distribution function 
of the emitting particles, rendering it difficult to distinguish between them on the basis 
of observations. However, the geometry of the magnetic field, which is the crucial factor 
determining the polarization properties, is constrained in the striped model to be close to 
that of the analytic asymptotic solution presented by Bogovalov (1999). We have therefore 
presented detailed computations of the polarization properties of the pulses expected in 
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this scenario. These possess the characteristic property, unique amongst currently discussed 
models, that the electric vector of the off-pulse emission is aligned with the projection of 
the pulsar's rotation axis on the plane of the sky. This is in striking agreement with recent 
observations of the Crab pulsar. In addition the striped wind scenario naturally incorporates 
features of the phase-dependent properties of the polarization angle, degree of polarization 
and intensity that are also seen in the data. This underlines the need to develop the model 
further, in order to confront high-energy observations of the Crab and other pulsars. In 
particular, the manner in which magnetic energy is released into particles in the current sheet 
remains poorly understood and the link between the asymptotic magnetic field structure and 
the pulsar magnetosphere is obscure. 

We thank Gottfried Kanbach for providing us with the OPTIMA data and for help- 
ful discussions. This work was supported by a grant from the G.I.F., the German-Israeli 
Foundation for Scientific Research and Development. 
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Models Observations 
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Fig. 1. — Light curve of intensity, degree of polarization and position angle of the pulsed 
synchrotron emission obtained by numerical integration of the set of equations (2), and 
measurements of these quantities for the Crab pulsar (Kellner 2002; Kanbach et al. 2003). 
Models with Lorentz factor T = 20 (solid line, red online) and 50 (dotted line, cyan online) 
are shown. The assumed particle energy distribution index was p = 2, and the inclination of 
the line of sight equals the obliquity: a = £ = 60°. The position angle of the projection of 
the pulsar's rotation axis was set to 124° (Ng & Romani 2004). The bottom panels show the 
dependence on phase (= Q*r/(2irv) in Eq. (1)) of the assumed magnetic field components 
and the particle density in the comoving frame. The maximum values of B v and K are 
normalized to unity. 



